Saturation of nuclear matter and short-range correlations 
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A fully self-consistent treatment of short-range correlations in nuclear matter is presented. Differ- 
ent implementations of the determination of the nucleon spectral functions for different interactions 
are shown to be consistent with each other. The resulting saturation densities are closer to the em- 
pirical result when compared with (continuous choice) Brueckner-Hartree-Fock values. Arguments 
for the dominance of short-range correlations in determining the nuclear-matter saturation density 
are presented. A further survey of the role of long-range correlations suggests that the inclusion 
of pionic contributions to ring diagrams in nuclear matter leads to higher saturation densities than 
empirically observed. A possible resolution of the nuclear-matter saturation problem is suggested. 

PACS numbers: 21.65.+f, 21.30.Fe 



A correct description of the saturation properties of 
nuclear matter has remained an unresolved issue for a 
very long time. The Brueckner-Bethe-Goldstone (BBG) 
expansion supplies a converged result for the energy 
per particle in the relevant density range, for a given re- 
alistic interaction, at the level of three hole-line contribu- 
tions Such calculations fail to reproduce the empir- 
ical saturation properties which require a minimum in the 
equation of state at a density corresponding to a Fermi 
momentum, hp, of about 1.33 fm _1 with a binding en- 
ergy of about 16 MeV. The authors of Ref. obtain for 
the Argonne V14 interaction a saturation density cor- 
responding to 1.565 fm -1 with about the correct amount 
of binding. This corresponds to an overestimation of the 
empirical density by about 60% but appears completely 
consistent with corresponding variational calculations 5] 
for the same interaction. 

Several different remedies for this serious problem have 
been proposed over the years. The intrinsic structure 
of the nucleon and its related strong coupling to the A- 
isobar inevitably requires the consideration of three-body 
(or more-body) forces. When three-body forces are con- 
sidered in variational calculations it is possible to achieve 
better saturation properties only when an adhoc repul- 
sive short-range component of this three-body force is 
added |6|,[7j. ^ nas a ^ so been suggested that a relativistic 
treatment of the nucleon in the medium using a Dirac- 
Brueckner approach provides the necessary ingredients 
for a better description of saturation [a, IfJ, [lfj, |ll| . 

All many-body methods developed for nuclear mat- 
ter have focused on a proper treatment of short-range 
correlations (SRC) without the benefit of experimental 
information on the influence of these correlations on the 
properties of the nucleon in the medium. This influence 
can now be clearly identified by considering recent re- 
sults from (e,e'p) reactions [TJ, LLil LLJ] and theoretical 
calculations of the nucleon spectral function in nuclear 
matter [lil Hfl Il7j. A recent analysis of the (e,e'p) re- 



action on 208 Pb in a wide range of missing energies and 
for missing momenta below 270 MeV/c yields informa- 
tion on the occupation numbers of all the deeply-bound 
proton orbitals. These data indicate that all these or- 
bitals are depleted by the same amount of about 15% . 
These occupation numbers are associated with the orbits 
which yield an accurate fit to the (e,e'p) cross section. 
The properties of these occupation numbers suggest that 
the main effect of the global depletion of these mean- 
field orbitals is due to SRC. Indeed, the effect of the 
coupling of hole states to low-lying collective excitations 
only affects occupation numbers of states in the immedi- 
ate vicinity of the Fermi energy |l9| . In addition, nuclear 
matter momentum distributions display such an overall 
global de plet ion due to short-range and tensor correla- 
tions |l7J, |2fJ, |2Jj . The latter results formed the basis of 
the now corroborated prediction [52, for the occupa- 
tion numbers in 208 Pb (lif. 

Most of this depleted single-particle (sp) strength is 
located at energies more than 100 MeV above the Fermi 
energy [rH l20t 122^ . This appearance of strength at high 
energy is another important aspect of the influence of 
short-range and tensor correlations. Yet another char- 
acteristic feature of these SRC is that this depletion of 
the sp strength must be compensated by the admix- 
ture of a corresponding number of particles with high- 
momentum components. These high-momentum compo- 
nents have not yet been unambiguously identified but 
are currently studied experimentally |24j . Solid theoreti- 
cal arguments |2?3 j and calculations clearly pinpoint this 
strength at high excitation energy in the hole spectrum 
both for nuclear matter |l7j and finite nuclei |26j. In- 
deed, experiment confirms that no substantial admixture 
of these high-momentum components is observed in the 
vicinity of the Fermi energy }27j . 

We now present an argument showing that SRC are 
the dominant factor in determining the empirical satu- 
ration density of nuclear matter. We recall that elastic 
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electron scattering from 208 Pb [23 accurately determines 
the value of the central charge density in this nucleus. 
By multiplying this number by A/Z one obtains the rel- 
evant central density of heavy nuclei, corresponding to 
0.16 nucleons/fm 3 or Uf — 1-33 fm _1 . Since the presence 
of nucleons at the center of a heavy nucleus is confined to 
s-wave nucleons, and, as discussed above, their depletion 
is dominated by SRC, one may therefore conclude that 
the same is true for the actual value of the empirical satu- 
ration density of nuclear matter. While this argument is 
particularly appropriate for the deeply bound lsj/2 and 
2s]/2 protons, it continues to hold to a large extent for 
the 3s 1/2 protons which are depleted predominantly by 
short-range effects (up to 15%) and by at most 10% due 
to long-range correlations [IjS, |23| . These considerations 
demonstrate clearly that one may expect SRC to have 
a decisive influence on the actual value of the nuclear- 
matter saturation density. 

High-momentum components due to SRC also have 
a considerable impact on the binding energy of nuclear 
matter. This result can be inferred from the energy sum 
rule [13 
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where p = is the density. Eq. JTJ illustrates the 

link between the energy of the system and the hole spec- 
tral function, Sh{k,oo). Results for the momentum dis- 
tribution and true potential energy based on the spectral 
function show that enhancements as large as 200% for 
the kinetic and potential energy over the mean-field val- 
ues can be obtained for both nuclear matter |2(J and 
finite nuclei [2(|. These large attractive contributions to 
the potential energy of nuclear matter are mainly from 
weighing the high-momentum components in the spec- 
tral function with large negative energies in Eq. • The 
location of these high-momentum components as a func- 
tion of energy is therefore an important ingredient in the 
determination of the energy per particle as a function of 
density. So far, the determination of this location has 
relied only on quasiparticle properties in the construc- 
tion of the self-energy. A self-consistent determination of 
the spectral function including the location of these high- 
momentum components therefore includes the dominant 
physics of SRC in the description of nuclear matter and 
is consistent with the experimental observations of the 
nucleoli spectral function in nuclei. 

Such a determination requires the solution of the lad- 
der equation for the effective interaction in the medium 

( q \r J J T (K,n)W) = ( q \v/ e ? T W) (2) 
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where a notation with relative momenta p, q, q' and the 
conserved total momentum K has been used. The prop- 
agator gj 1 in Eq. J5J has been obtained by an angle- 
averaging procedure of the noninteracting two-particle 
propagator 
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in order to allow a partial wave decomposition of the 
ladder equation. Note that sp momenta k, k' are used 
in Eq. ©. In turn, the spectral functions are needed to 
determine Eq. © and thereby the effective interaction T 
through Eq. J2J. They can be obtained from the imag- 
inary part of the sp propagator which solves the Dyson 
equation 

g(fc, to) = g<°> (k, to) + g<°> (fc, to)E(k, to)g(k, to), (4) 

where the self-energy £ includes the contribution of SRC 
through T, to complete the self-consistency loop. 

The implementation of this self-consistency scheme is 
numerically quite involved and has been attempted by 
several groups 0, 113, HU |3 HE EH- In the present 
paper two different approaches have been used to gen- 
erate results for different interactions. In the continuous 
scheme, a representation of the imaginary part of the self- 
energy in terms of four gaussians is used to completely de- 
scribe the sp propagator. The parameters of these gaus- 
sians are then determined self-consistently for the 
Reid potential |33 by solving Eq. J5J with the convolu- 
tion of spectral functions in Eq. © as input, constructing 
the self-energy, and then solving the Dyson equation 10}. 
In the discrete scheme we used a representation of the 
propagator in terms of three discrete poles |3^|, which 
avoids a full continuum solution of Eq. J5J. The latter 
approach is equivalent to a continuous version as far as 
the energy per particle is concerned, since it requires a 
reproduction of the relevant energy-weighted moments 
of the hole and particle spectral function [3^]. This is 
substantiated by comparing the results of this discrete 
scheme with the results of the continuous self-consistency 
scheme used in Ref. [34| for the Mongan-type separable 
interaction |33 and in Ref. [3^ for the separable Paris 
interaction 40]. We find that the binding energies cor- 
respond to within 5% over the relevant kp range around 
the minimum, and moreover that the location of the min- 
imum agrees to within 3%. 

In Fig^ we report the saturation points obtained 
within the discrete scheme of Ref. |33J for the updated 
Reid potential (Reid93), the Nijml and Nijmll inter- 
action [3^ and the separable Paris interaction [4o| . 
The results demonstrate an important and systematic 
change of the saturation properties with respect to con- 
tinuous choice Brueckner-Hartree-Fock (ccBHF) calcula- 
tions, leading to about 4-6 MeV less binding, and reduced 
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FIG. 1: Nuclear matter saturation points calculated with 
various realistic NN interactions. The open symbols refer to 
continuous choice Brueckner-Hartree-Fock results. The filled 
symbols refer to self-consistent results and represent satura- 
tion points calculated in the discrete scheme, except for the 
Reid68 interaction where the binding energy at two densities 
is shown in the continuous scheme. 



values of the saturation density, closer to the empirical 
one. Such a trend is entirely consistent with the obser- 
vations in Ref. 36] made for a separable NN interaction, 
and is now extended to more realistic (non-separable) 
interactions. 

The discrete scheme could not be used with the 
original Reid (Rcid68) potential because of its slow de- 
cay in momentum space, but some results are available in 
the continuous scheme of Ref. • In Fig^ the binding 
energy is shown at two densities (kF = 1-33 and 1.45 
fm _1 ); the error bars are an estimate of the remain- 
ing uncertainty due to incomplete convergence and the 
non-selfconsistent treatment of some higher order partial 
waves :32]. The results again seem to indicate a substan- 
tial shift in the saturation density for the Reid68 poten- 
tial, from the ccBHF value of about 1.6 fm _1 , to a value 
below 1.45 fm _1 , without seriously underbinding nuclear 
matter. 

The present self-consistent treatment of SRC (sc- 
SRC) differs in two main aspects from the ccBHF ap- 
proach. Firstly, hole and particle lines are treated 
on an equal footing, thereby ensuring thermodynamic 
consistency 34J . Intermediate hole-hole propagation in 
the ladder diagrams is included to all orders. This feature 
provides, compared to ccBHF, a substantial repulsive ef- 
fect in the k < kF contribution to Eq. JQl, and comes 
primarily from an upward shift of the quasi-particle en- 
ergy spectrum as a result of including uj < tp contribu- 
tions to the imaginary part of the self-energy. The ef- 



fect increases with density, and is the dominant factor in 
the observed shift of the saturation point. Secondly, the 
realistic spectral functions, generated through Eqs.(j2EJ 
and used in the evaluation of the in-medium interaction 
r and self-energy S, are in agreement with experimen- 
tal information obtained from (e,e'p) reactions. For the 
Reid93 interaction at kF = 1.37 fm" 1 we find z = 0.74 
for the quasiparticle strength at the Fermi momentum, 
whereas the hole strength for p = 0, integrated up to 
100 MeV missing energy, equals 83%; similar values are 
found for the other interactions. The depletion of the 
quasiparticle peaks is primarily important to suppress 
unrealistically large pairing instabilities around normal 
density. The improved treatment of the high-momentum 
components does affect the binding energy, through the 
k > kF contribution to Eq. JQ). This feature, studied in 
|33|. provides a sizeable attraction, but is smaller than 
the afore- mentioned repulsive effect. 

The inclusion of /i/i-propagation in scSRC also leads 
to a somewhat stiffer equation-of-state than in ccBHF. A 
recent analysis of the giant monopole resonance in heavy 
nuclei |4l| yields an experimental estimate K nm — 210 ± 
30MeV for the nuclear matter compression modulus, 
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At the saturation points in FigQ] we find ccBHF values 



Kr, 



154 MeV for Reid93 and K„ 



148 MeV for 



the separable Paris interaction, which are enhanced to 



Kr, 



177 MeV and K n 



216 MeV, respectively, in 



our scSRC calculation. These values agree reasonably 
well with the experimental estimate. Note that reason- 
able values for K imply that the Reid68 energies in FigQ] 
may still deviate by 1-1.5 MeV from numerically exact 
scSRC values, as indicated by the error bars. 

The present results indicate that a sophisticated treat- 
ment of SRC lowers the ccBHF saturation densities, 
bringing them closer to the empirical one. It remains to 
be understood why apparently converged hole-line calcu- 
lations Q yield higher saturation densities. The three 
hole-line terms obtained in Ref. |3| indicate reasonable 
convergence properties compared to the two hole-line 
contribution. One may therefore assume that these re- 
sults provide an accurate representation of the energy per 
particle of nuclear matter as a function of density for the 
case of nonrelativistic nucleons and two-body forces. At 
this point it is useful to identify an underlying assump- 
tion when the nuclear-matter problem is posed 32] . This 
assumption asserts that the influence of long-range cor- 
relations in finite nuclei and nuclear matter are commen- 
surate. We'd like to point out that this underlying as- 
sumption is questionable. Three hole-line contributions 
include a third-order ring diagram characteristic of long- 
range correlations. The effect of long-range correlations 
on nuclear saturation properties is sizeable, as shown by 
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the results for three- and four-body ring diagrams cal- 
culated in Ref. H3 (see also Refs. 0, Q). The results 
of Ref. ^3 demonstrate that such ring-diagram terms 
are dominated by attractive contributions involving pion 
quantum numbers propagating around the rings, and in- 
crease in importance with increasing density. Such long- 
range pion-exchange contributions to the binding energy 
appear due to the possibility to coherently sample the at- 
tractive interaction in a given ring diagram at momenta 
q above 0.7 fm^ 1 . This feature is related to momentum 
conservation in nuclear matter and is not available in fi- 
nite nuclei, in which no such collective pion-degrees of 
freedom are actually observed |4j|. It seems therefore 
reasonable to call into question the relevance of these 
coherent long-range pion-exchange contributions to the 
binding energy per particle since their behavior is so 
markedly different in finite and infinite systems. One 
may consider the salient difference of the ratio of spin- 
longitudinal and spin-transverse response function in nu- 
clear matter and finite nuclei as another indication of 
the relevance of our suggestion 0] . We also like to point 
out that experimental information of these response func- 
tions |45|, llfl, |42| suggests no characteristic enhancement 
of the (pionic) spin-longitudinal response as expected on 
the basis of nuclear matter calculations. 

Clearly, the assertion that long-range pion-exchange 
contributions to the energy per particle need not be con- 
sidered in explaining nuclear saturation properties, needs 
to be further investigated. At this point it appears that 
a fully self-consistent treatment of SRC has substantially 
different saturation properties than a conventional (con- 
tinuous choice) Brueckner-Hartree-Fock treatment, and 
is capable of yielding saturation densities close to the 
empirical one. 
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